Coupling of frustrated Ising spins to magnetic cycloid in multiferroic TbMnOs 
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We report on difTraction measurements on multiferroic TbMnOa which demonstrate that the Tb- 
and Mn-magnetic orders are coupled below the ferroelectric transition Tfe ~ 28 K. For T < Tfe 
the magnetic propagation vectors (r) for Tb and Mn are locked so that t'^^— t"^", while below 
T'n'^= 7 K we find that t'^^ and t*^" lock-in to rational values of 3/7 b* and 2/7 b*, respectively, 
and obey the relation 3t"^^- r^"= 1. We explain this novel matching of wave vectors within the 
frustrated ANNNI model coupled to a periodic external field produced by the Mn-spin order. The 
T^^= T^" behavior is recovered when Tb magnetization is small, while the t'^^= 3/7 regime is 
stabilized at low temperatures by a peculiar arrangement of domain walls in the ordered state of 
Ising-like Tb spins. 

PACS numbers: 61.12.Ld, 61.10.-i, 75.30.Kz, 75.47.Lx, 75.80.+q 



While ferroelectricity (FE) and magnetism are chemi- 
cally incompatible [l[ , it has recently been shown that in- 
version and time-reversal symmetry can be broken simul- 
taneously if magnetic spins order in a cycloidal arrange- 
ment @,y,l3|- This mechanism is realized in i?Mn03 per- 
ovskite manganites (i? =Tb and Dy) where the magnetic 
ordering of Mn-spins forms a cycloid within the 6c— plane. 
Accordingly, the spontaneous ferroelectric polarization, 



Ps, is given by Ps oc myUi^e^ x r), where my and respectively, and hold the relationship St 



are the components of the Mn spins in the cycloid along 
the b— and c— axis respectively, ex is a unit vector along 
the axis of rotation of Mn-spins and r is the magnetic 
propagation vector 0] . Recent studies showed that direct 
and indirect contributions of rare earth spins to electric 
polarization can also be very important. In DyMnOa, for 
example, the tripling of Pg is quantitatively accounted 
by the effective enhancement of my by the ordering of 
Dy-spins along the 6— axis 6]. Further influence of R- 
ions on the properties of multiferroics are found in the 
flop of the cycloid into the a6— plane (and consequently, 
Ps||a) for the case of the non-magnetic R =(Eu,Y) [7[. 
Although the cycloidal ordering of Mn-spins drives mul- 
tiferroicity, i?-ions strongly modulate it and thus signifi- 
cantly influence multiferroic properties. 

Similar effects are also evident in the i?Mn205 class of 
multiferroics despite the fact that here multiferroic be- 
havior may arise from an exchange striction mechanism 
Q. For example, in i? = Tb, increases in the low 
temperature incommensurate phase, while it decreases 
and reaches a negative value for non-magnetic R = Y 
lol m . This difference in behavior results from the in- 
fluence of Tb-spin ordering upon the magnetic ordering 
of Mn-spins |8l|. It is therefore clear that irrespective 
of the mechanism that drives multiferroic behavior, the 
magnetic coupling between R- and Mn-spins {Jmu-r) 
needs to be understood in order to arrive at a detailed 



and quantitative model of multiferroics and understand 
ways to improve their performance. 

In this Letter we report on diffraction measurements 
which demonstrate that the Tb- and Mn- magnetic order- 
ing in multiferroic TbMnOs remain coupled for T < Tfe. 
While below TpE it is known that Tb-spins are induced 
to order with t^^= t'^" 0, below 7 K we find that x^^ 
t'^" lock-in to rational values of 3/7 b* and 2/7 b*, 

1. We 



and 



Tb 



explain the novel matching of Tb and Mn wave vectors 
within the frustrated ANNNI model coupled to a peri- 
odic external field produced by the Mn-spin ordering. 
The t'^^— r'^" behavior is recovered and we show that 
surprisingly the t'^^= 3/7 regime is stabilized by an opti- 
mal Tb spin-density wave ordering with 6 domain walls, 
superimposed on the t'^"= 2/7 Mn-ordering. Our phase 
diagram not only reproduces the experimentally observed 
value of t"^^ but also explains magnetic orde ring for R- 
spins that is realized for R = Dy and Ho [l^, [TsTIT^I . 

High quality single crystals of TbMnOs were obtained 
using the floating zone technique and characterized us- 
ing magnetization measurements with excellent agree- 
ment with pubhshed data 15, [3]. Single crystal neu- 
tron diffraction measurements were performed on two 
TbMnOa crystals using A = 2.43 A, at the El triple-axis 
spectrometer operated in 2-axis mode and located at the 
BENSC facility of the Hahn-Meitner-Institute (HMI). To 
directly measure the magnetic ordering of Tb-spins and 
compare it to the ordering of Mn-spins obtained from 
neutron diffraction we have used X-ray resonant mag- 
netic scattering (XRMS) at the Tb L2,3 edges and 12.4 
keV non-resonant X-ray diffraction (XRD) at the 7 T 
multipole wiggler beamline MAGS, operated by the HMI 
at the synchrotron source BESSY [l3|. For the XRMS 
measurements, linear polarization analysis was carried 
out using HOPG(006) and Al(lll) crystals, respectively. 
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low T]^" 



The a— axis was kept parallel to the scattering plane. 
Additional XRMS measurements with applied magnetic 
field H\\a and a perpendicular to the scattering plane 
were performed at the XmaS beamline at the ESRF, 
Grenoble. 

Our single crystal neutron diffraction experiments be- 
largely confirm previous measurements [3, fisl . 
Il9j . We find a set of magnetic satellites in Brillouin zones 
(hkl) with extinction conditions h + k = even, I = odd 
that arise exclusively from the ordering of Mn-spins. Us- 
ing the convention of Bertaut these are known as A-type 
reflections The magnetic ordering of Mn-spins is 

incommensurate and characterized by the propagation 
vector T^^" ~ 0.285 b*. For T < Tfe = 28 K, i.e., 
in the ferroelectic phase, reflections in all G, C, F-type 
Brillouin zones are observed 



25| 



These reflection arise 
from the polarization of the Tb-spins along the a— axis 
U,^. Below T^b = 7K, a strong increase of intensities 
of G-, F- and C-type satellites with t^" = 0.277(2) b* 
is observed and much more intense magnetic satellites 
appear at positions characterized by a propagation vec- 
tor t'^^ = 0.426(2) b* and its odd harmonics. The latter 
wave vector has been interpreted as to describe a separate 
ordering of Tb-spins below T^^. However, the tempera- 
ture dependence of r'^" reflections suggest that Tb-spins 
may also continue to order in this regime with r^". 

In contrast to previous investigations we separate the 
contribution of Tb spin ordering from that of Mn to the 
various superlattice reflections by employing XRMS. In 
Fig. [T^,b we show contour plots of the temperature de- 
pendences of selected reflections. Below Tfe, a number 
of relatively weak resonant (0, fc ± r'^",/) reflections of 
G- and C-type have been observed in the cr — tt' polariza- 
tion channel, representing the magnetic ordering of Tb 
with the same propagation vector as the Mn ordering, 
^Tb ^ ^Mn ^ 0.274(1) at 27 K. These reflections are 
observed below TpE and their intensity increases with 
concave curvature on cooling - a behavior typical for an 
induced moment. With decreasing temperature below 
T^^ the value of t^^ shifts to a higher value 0.286(1) 
but more surprisingly we find no decrease in its inten- 
sity in the magnetic channel, a behavior also noted in 
our neutron diffraction data. Rather we observe a con- 
tinuous increase suggesting that the magnitude of Tb- 
spins at this wave vector increases below T^^. This ef- 
fect is accompanied by the appearance of much stronger 
(0, k ± tJ^, I) reflections, with A- and F-type being dom- 
inant, corresponding to a separate ordering of Tb with 
tJ^ = 0.4285(10). The observed sequence of magnetic 
transitions corresponds to Tb magnetic order with t^^ — 
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FIG. 1: (Color online) Temperature dependence of charac- 
teristic magnetic and lattice distortion reflections in the form 
of a color intensity plot measured using XRMS and XRD. In 
panels (a) and (b) we show cr — tt' data measured from the 
t'^^ and r*^" peaks in the (0 k 2) zone with the X-ray energy 
tuned to the Tb L2 edge. In panels (c) and (d) we shown 
XRD data of the r"^" and 2r^" reflections measured in the 
(0 k 3) and (0 k 2) Brillouin zones, respectively. The r^" 
reflections in panel (c) is 2 x lO'' smaller in intensity than 
the neighboring (023) reflections and represent only magnetic 
scattering. Below T^^ the weak (0,2 -I- r*^", 3) reflection is 
superposed by a 100 times stronger reflection from a strucural 
distortion, —4t'^^, related to the Tb ordering. The intensity 
scale in panel (c) is broken close to T^^ and a linear back- 
ground has been subtracted. The 2r^" in panel (d) is the 
expected second harmonic structural reflection 



T < Tfe), we employed XRD. In Fig.[TJ;,d we show tem- 
perature dependent measurements of the non-resonant 
A-type reflection (0, 2 + t^", 3) (which has only 2 x 10"^ 
the intensity of the neighboring (0 2 3) reflection) and the 
second harmonic reflection (0, 2 + 2t'^", 2), representing 
directly the Mn magnetic ordering and the accompanied 
lattice distortion, respectively. Below T^^, both peaks 
shift simultaneously showing that indeed t'^" shifts from 
0.274 at T^^ to 0.286 just below it. Having estabhshed 
the behavior of the Tb-ordering, from here on we simplify 
our notation by taking t'^^{= t]'^) and r'^"(= t'[^). 

We note that below T^^ our x-ray scattering data re- 
veals that to an accuracy of 10^'' the two wave vectors 
are locked-in to the rational values t'^^ — 3/7 b* and 
_Mn _ 2/7 b*. These values are slightly larger than 



0.274 for T^^< T < Tfe and Tb magnetic or- those reported by neutron diffraction @, where the 



0.286 below T^^. 



der with essentially two wave vectors rj*^ =0.4285 and shift of r'^" across T^^ is much smaller than observed 

here [2^. The lock-in of t^^ and r^" to values of ratio- 
nal fractions reveals that below T^^ a portion of Tb-spins 
does not in fact order with yet a separate wave vector, a 
behavior that would be unusual. Rather, the two wave 



To confirm that the shift of t]^'"=0.274: 0.286 at 
(on cooling) also reflects the behavior of the Mn 



rpTh 

sublattice (i.e. the coupling t^^ = r^" is valid for all 
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FIG. 2: (color online) Average values of R spins (red circles) 
described by the frustrated Ising model in periodically varying 
external field (dotted blue line) at zero temperature (a) and 
at high temperatures (b). In the latter case r^'' — r^" = 2/7, 
while at T = 0, t'^'° = 3/7. 
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FIG. 3: (color online) Zero temperature phase diagram of 
model lU. The phases are numbered by the wave vector r'^^ 
of the main spin modulation: r^'' = 2/7 (yellow), r^'' = 3/7 
(red), r'^'" = 1/2 (green). 



vectors are related by Sr"""*^ - r'^" — 1, implying a single 
phase Tb magnetic structure below T'P^. This relation 
does not only hold for the present x-ray scattering data 
but also for our neutron diffraction measurements, and 
others reported 0, [l^l, to an accuracy of better than 
0.002. Therefore, this relation appears to hold irrespec- 
tive if the values of r^" and t'^^ are exactly 2/7 and 
3/7. The behavior we describe here clearly points to an 
intimate relation between the Mn and Tb magnetic or- 
derings mediated by JMn-Tb exchange and suggests that 
the ordering of Tb spins as well as the ground state wave 
vectors are determined by this interaction. 

The 3t"^^ - r^" = 1 ordering can be modulated by ap- 
plying a magnetic field along the a— axis. Using XRMS 
measurements with H\\a at 2 K (not shown), we find 
a first order phase transition at 1 T characterized by 
vanishing of t'^^ = 3/7 reflections and the appearance 
of new reflections at t"^^ = 0.361(1) b*, and a shift of 
^Mn 0.285(1) down to 0.276(1). This behavior sug- 
gests that now the two wave vectors become related by 

Phenomenologically, the zero fleld behavior results 
from the quartic coupling gL^^^LMn between the Tb and 
Mn magnetic order parameters, which favors a relation 
of the type St"^^ ± t^" = n (n is an integer number) be- 
tween the wave vectors of the two spin-density waves and. 



Tb 



1 



in particular, gives the observed relation 3r 
for TbMnOa at iJ = 0. In applied magnetic field H, the 
Tb spin-density wave acquires a homogeneous component 
xH and the quartic term effectively leads to the coupling 
3gxH L^^Lun favoring 2r'^'^ ± r'^" = n, explaining the 
switch of the relation between the Tb- and Mn- spins 
observed for H\\a. 

This phenomenological description applies close to the 



Tb ordering temperature. At zero temperature Tb-spins 
are Ising-like (22,, 23] and cannot be described by a small 
number of Fourier harmonics. To understand the inter- 
play between the R and Mn ordcrings we consider the 
one-dimensional ANNNI Ising model. 



E ■ 



E 

y 



[Jaya. 



y+l/2 



J'ayay + l — \Uy COS^Ql/ + I 



(1) 

where the Ising variables Cj, = ±1 describe R spins at 
the sites y = 0, 1/2, 1, 3/2, 2, . . ., coupled to the period- 
ically varying 'external field' induced by Mn spins, so 
that Q = r^" and A oc Jivin-Tb- The unit cell of the 
chain contains two spins (one with an integer and one 
with an half-integer coordinate) corresponding to two R 
ions in each a6— layer of the unit cell of i?Mn03. For fer- 
romagnetic nearest- neighbor coupling, J < 0, antiferro- 
magnetic next-nearest-coupling, J' > 0, and A = we re- 
cover the usual ANNNI model, which has the FM ground 



state (r 



Tb 



0) for K 



\J\ 



< 



and the |f||-state 



{t^^ = 1/2), for k > i. At the frustration point k = 1/2 
the FM domain wall (DW) energy Edw = -2(2J' + J) 
is zero and the ground state configuration can have any 
number of domain walls provided that they do not oc- 
cupy nearest-neighbor sites 24 1 . 

We now compare energies of the states with t^^ = 2/7, 
3/7 and 1/2 corresponding to, respectively, 4, 6, and 7 
DWs in the interval of 7 unit cells (14 spins). For k > 1/2 
the DW energy is negative and the ttii-state having the 
highest number of DWs has the lowest energy, if only 
interactions between Tb-spins are considered. The en- 
ergy of the Tb-Mn interaction in "f m-statc is zero, while 
for the state with t^'^ = r'^''" = 2/7 it is -0.642A per 
spin. Surprisingly, the state with t^^ — 3/7 also gains a 
substantial energy due to the interaction with Mn spins. 
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which for the best possible arrangement of 6 DWs, shown 
in Fig. [2^, is — 0.218A per spin. 

The resulting zero temperature phase diagram is shown 
in Fig. O For the weak R-Mn coupling regime, the TTii" 
state has the lowest energy and the wave vectors of R 
and Mn orderings are decoupled as is the case for R 
= Dy [13] ■ For large A, = down to zero tem- 
perature as found for R = Ho 13|, [ij]. In the narrow 
sector 0.656 < jj^+j < 0.674, Ising spins order in the 
ground state with t^^ = 3/7, as observed in TbMnOa. 
This is a rather unexpected result since the 3/7-state 
is stabilized by the modulation of Mn-spins with a dif- 
ferent wave vector. Numerical mean field studies of the 
three-dimensional version of model ([T]) for parameters fa- 
voring the 3/7-state show that at high temperatures the 
wave vector of the Tb modulation equals that of Mn- 
spins (see Fig. [^b), while below T^^ the intensity of the 
t"^^ = 3/7 peak is several times higher than that of the 
t"^^ = 2/7 peak, in agreement with our diffraction data. 
The small extent of the 3/7-region in the phase diagram 
of our model suggests that other mechanisms, e.g. longer 
range interactions between i?-spins and the deformation 
of the Mn spin-density wave in response to the ordering 
Tb-spins, may contribute to stabilization of this phase. 

Our results uncovered a strong coupling between Tb- 
and Mn-spin sublattices in multiferroic TbMnOs that is 
allowed because of the FE nature of this compound at 
low temperatures. Indeed, Mn magnetic ordering alter- 
nating in the c-direction should produce zero exchange 
fields at the i?-ion as noted for HoMnOa in Ref. 13| . We 



believe that this argument is true above the FE transition 
in i? = Tb, Dy and Ho, where consequently no induced 
moment is found. In turn, the observation of induced 
ordering below the FE transition can be considered as 
strong evidence of symmetry breaking that would allow 
i?-moments to be induced as the spontaneous electric po- 
larization occurs. 

In summary, our work demonstrates that the Jmii-r 
interaction is a critical ingredient in our understand- 
ing of the magnetism and physical properties of mul- 
tiferroic manganites with magnetic i?-ions, highlighting 
the relevance of our model. In the temperature range 
< T < TpE the coupling of R- and Mn-spin order is 
reflected in = r'^", a behavior that is now confirmed 
for R =Tb, Dy and Ho [3, Hill- For T <T^, R=Ro 

Ho 

Our model 



represent a strong coupling regime as r'^"= r"" is pre- 



served to the lowest temperatures [1 
reproduces the behavior for i?=Dy for a weak coupling 
regime as we recover the commensurate t^'^=1/2 order- 
ing [1, [13]. However, the case of i? = Tb is more inter- 
esting as it represents an intermediate coupling regime 
with complex behavior. In particular, below T^^^, when 
Jxb-Tb becomes important, we find that the Tb- and 
Mn- orderings remain coupled through the matching of 
their wave vectors, St"^^ — t^" = 1. The adjustment of 
Ising-like Tb spins to periodic Mn ordering results in a 



rather complex shape of the Tb spin-density wave at low 
temperatures. 

Finally, the coupling between R- and Mn-sublattices 
in .RMnOa with magnetic R = Tb, Dy and Ho may have 
some implications for magnetically-induced ferroelectric- 
ity in these multiferroics. Referring to Eui-^^YajMnOa 
with nonmagnetic R and Mn-cycloid within the a5-plane 
[tI, the stabilization of the 6c-cycloid for R = Tb and 
Dy may be caused by the energy gain due to the Tb-Mn 
exchange in the FE state, which allows non-zero Mn ex- 
change fields at R sites. The counterintuitive magnetic 
field behavior for R = Tb and Dy, e.g. the cycloid fiop 
from the be- to a6-plane in H\\a, can then result from 
reduction of the R-Mn interaction when i?-spins become 
aligned ferromagnetically by the magnetic field. 
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Laurence Bouchenoire from ESRF for assistance during 
the experimental work. M.M. thanks HMI for hospitality. 



[9: 

[lo; 

[11 

[12; 

[is: 

[14: 

[is: 
[16: 



N.A. Hill, J.Phys. Chem. B 104, 6694 (2000). 
M. Kenzelmann, A. B. Harris, S. Jonas, C. Broholm, 
J. Schefer, S. B. Kim, C. L. Zhang, S.-W. Cheong, O. 
P. Vajk, J. W. Lynn, Phys. Rev. Lett. 95, 087206 (2005). 
M. Mostovoy, Phys. Rev. Lett. 96, 067601 (2006). 

H. Katsura, N. Nagaosa, A. V. Balatsky, Phys. Rev. Lett. 
95, 057205 (2005). 

T. Goto, T. Kimura, G. Lawes, A. P. Ramirez, Y. Tokura, 
Phys. Rev. Lett. 92, 257201 (2004). 

O. Prokhnenko, R. Feyerherm, E. Dudzik, S. Landsgesell 
N. Aliouane, L. C. Chapon, D. N. Argyriou, Phys. Rev. 
Lett. 98, 057206 (2007). 

J. Hemberger, F. Schrettle, A. Pimenov, P. Lunken- 
heimer, V. Yu. Ivanov, A. A. Mukhin, A. M. Balbashov, 
A. LoidI, Phys. Rev. B 75, 035118 (2007). 
L.C. Chapon, G. R. Blake, M. J. Gutmann, S. Park, 
N. Hur, P. G. Radaelli, S.-W. Cheong, Phys. Rev. Lett. 
93, 177402 (2004). 

L.C. Chapon, P. G. RadaeUi, G. R. Blake, S. Park, S.- 
W. Cheong, Phys. Rev. Lett. 96, 097601 (2006). 

I. Kagomiya, S. Matsumoto, K. Kohn, Y. Fukuda, 
T. Shoubu, H. Kimura, Y. Noda, N. Ikeda, Ferroelectrics 
286, 889 (2003). 

N. Hur, S. Park, P. A. Sharma, J. Ahn, S. Guha, S.- 

W. Cheong, Nature 429, 392 (2004). 

R. Feyerherm, E. Dudzik, N. Aliouane, D. N. Argyriou, 

Phys. Rev. B 73, 180401 (R) (2006). 

A. Munoz, M. T. Casais, J. A. Alonso, M. J. Martmez- 

Lope, J. L. Martinez, M. T. Fernandez-Di'az, Inorg. 

Chem. 40, 1020 (2001). 

H.W. Brinks, J. Rodn'guez-Carvajal, H. Fjellvag, 
A. Kjekshus, B.C. Hauback, Phys. Rev. B 63, 094411 
(2001). 

T. Kimura, T. Goto, H. Shintani, K. Ishizaka, T. Arima, 
Y. Tokura, Nature 426, 55 (2003). 

T. Kimura, G. Lawes, T. Goto, Y. Tokura, A. P. Ramirez, 
Phys. Rev. B 71, 224425 (2005). 



5 



[17] E. Dudzik, R. Fcyorhcrm, R. Signorato, C. Zilkens, J. 

Synchrotron Rad. 13, 421 (2006). 
[18] S. Quczcl, F. Tchcou, J. Rossat-Mignod, G. Quezel, 

E. Roudaut, Physica B&C 86-88, 916 (1977). 
[19] R. Kajimoto, H. Yoshizawa, H. Shintani, T. Kimura, 

Y. Tokura, Phys. Rev. B 70, 012401 (2004). 
[20] E.F. Bertaut, in Magnetism, edited by G. T. Rado and 

H. Suhl, (Academic, New York, 1963) III, 149 (1963). 
[21] N. Aliouanc, D. N. Argyriou, K. Habicht, S. Landsgesell, 

K. Rule, in preparation (2007). 
[22] H.J. Bielen, J. Mareschal, J. Sivardiere, Z. Angew. Phys. 

23, 243 (1967). 



[23] J.E. Bourcc, J. Hammann, J. Phys. Paris 36, 391 (1975). 

[24] M.E. Fisher, W. Seiko, Phys. Rev. Lett. 44, 1502 (1980). 

[25] G: h + k = odd, I = odd; C: h + k = odd, I = even; F: 
h + k = even, I = even. 

[26] This discrepancy appears to be related to surface strain 
effects, since the present x-ray scattering measurements 
probe a thin surface-near region of a few jim only and the 
neutron measurements can bo reproduced using 100 keV 
x-rays (J. Strempfer et al., manuscript in preparation). 
Our basic conclusions are unaffected by this. 



